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Abstract The functional renormalization group provides an efficient description of 
the interplay and competition of correlations on different energy scales in interact- 
ing Fermi systems. An exact hierarchy of flow equations yields the gradual evolu- 
tion from a microscopic model Hamiltonian to the effective action as a function of 
a continuously decreasing energy cutoff. Practical implementations rely on suitable 
truncations of the hierarchy, which capture nonuniversal properties at higher energy 
scales in addition to the universal low-energy asymptotics. As a specific example we 
study transport properties through a single-level quantum dot coupled to Fermi liq- 
uid leads. In particular, we focus on the temperature T = gate voltage dependence 
of the linear conductance. A comparison with exact results shows that the functional 
renormalization group approach captures the broad resonance plateau as well as the 
emergence of the Kondo scale. It can be easily extended to more complex setups of 
quantum dots. 



1 Introduction 

The coupling of a quantum dot with spin degenerate levels and local Coulomb in- 
teraction (modeled, e.g., by the single-impurity Anderson model; SIAM) to metallic 
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leads gives rise to Kondo physics [1]. At low temperatures and for sufficiently high 
barriers the local Coulomb repulsion U leads to a broad resonance plateau in the 
linear conductance G of such a setup as a function of a gate voltage V g which lin- 
early shifts the level positions [2, 3, 4, 5, 6]. It replaces the Lorentzian found for 
noninteracting electrons. On resonance the dot is half-filled implying a local spin- 
5 degree of freedom responsible for the Kondo effect [1]. For the SIAM the zero 
temperature conductance is proportional to the one-particle spectral weight of the 
dot at the chemical potential [7]. The appearance of the plateau in the conductance 
is due to the pinning of the Kondo resonance in the spectral function at the chem- 
ical potential for — U /2 < V g < U/2 (here V g = corresponds to the half-filled 
dot case) [1, 6]. Kondo physics in transport through quantum dots was confirmed 
experimentally [8, 9], and theoretically using the Bethe ansatz [2, 6] and the numer- 
ical renormalization-group (NRG) technique [10, 11]. However, both methods can 
hardly be used to study more complex setups. In particular, the extension of the NRG 
to more complex geometries beyond single-level quantum-dot systems [5, 6, 12] is 
restricted by the computational complexity which increases sharply with the number 
of interacting degrees of freedom. Alternative methods which allow for a systematic 
investigation are therefore required. We here propose the functional renormalization 
group (fRG) approach to study low-temperature transport properties through meso- 
scopic systems with local Coulomb correlations. 

A particular challenge in the description of quantum dots is their distinct behav- 
ior on different energy scales, and the appearance of collective phenomena at new 
energy scales not manifest in the underlying microscopic model. An example of this 
is the Kondo effect where the interplay of the localized electron spin on the dot and 
the spin of the lead electrons leads to an exponentially (in U) small scale Tg . This 
diversity of scales cannot be captured by straightforward perturbation theory. One 
tool to cope with such systems is the renormalization group: by treating different 
energy scales successively, one can often find an efficient description for each one. 
We will here give an introduction to one particular variant, the fermionic fRG [13], 
which is formulated in terms of an exact hierarchy of coupled flow equations for 
the vertex functions as the energy scale is lowered. The flow starts directly from 
the microscopic model, thus including nonuniversal effects on higher energy scales 
from the outset, in contrast to effective field theories capturing only the asymptotic 
behavior. As the cutoff scale is lowered, fluctuations at lower energy scales are suc- 
cessively included until one finally arrives at the full effective action (the generating 
functional of the one-particle irreducible vertex functions) from which all physical 
properties can be extracted. This allows to control infrared singularities and com- 
peting instabilities in an unbiased way. Truncations of the flow equation hierarchy 
and suitable parametrizations of the frequency and momentum dependence of the 
vertex functions lead to powerful new approximation schemes, which are justified 
for weak renormalized interactions. A comparison with exact results shows that the 
fRG is remarkably accurate even for sizeable interactions. 

The outline of this article is as follows. In Sec. 2 we introduce the fRG formalism 
and derive the hierarchy of flow equations. The implementation of the fRG tech- 
nique for a quantum dot modeled by a SIAM is described in Sec. 3. In the following 
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we present results for the linear conductance including a comparison with exact so- 
lutions and discuss the emergence of the Kondo scale. We conclude in Sec. 4 with a 
summary and outlook on further applications and extensions of the present work. 



2 The fRG technique 

To make this article self-contained, this section aims to give a very short introduction 
into the fRG formalism by deriving the functional flow equations. Readers not so 
much interested in the formal beauty might jump ahead to the following Sec. 3 
where these flow equations are applied to a specific physical model. An introduction 
to the many -body formalism used here can be found, e.g., in [14]; for the details 
of the derivation of the functional flow equation see, e.g., [15, 16] and references 
therein. Exact functional flow equations were derived for bosonic field theories in 
[17, 18, 19] and for fermionic fields in the one-particle irreducible scheme in [13]. 



2.1 Green and vertex functions 

We consider a system of interacting fermions described by Grassmann variables \y, 
xfr, and an action 

S[\i/,y} = (y,C- 1 \if)-V[y,y} (1) 

with bare propagator 

C{K) = JT^T (2) 

IE - <5fc 

where the index K = (e,k, a) collects the Matsubara frequency e and the quantum 
numbers of a suitable single-particle basis set, e.g., momentum k and spin projection 
(7, and = £k — pi denotes the energy relative to the chemical potential. In order to 
understand the general structure of the flow equation it is useful to include a Nambu 
(particle/hole) index into K: then each directed fermion line denotes the propagation 
of either a particle or a hole [13]. The inner product implies a summation over 
these indices: for our diagonal propagator, (\jf,C~ l \l/) —Y.kWkC^^k- V[\j/,\j/] is 
an arbitrary many-body interaction; we will see a specific example for this below in 
Eq. (17). Connected Green functions can be obtained from the generating functional 
[14] 

^[77,77] = -log(.yK [Dyye^Me-^-^] (3) 



by taking derivatives with respect to the source fields 77 : 
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G m {K[,... 7 K' m ;Ki , . . . ,K m ) — - ■ ■ ■ Vkj, ¥k„, ■ ■ ■ Wi )conn 

gm gm 



8rj Kl ...8ri Km <?%/_. . .8r) K , 



n=n=o 



Equivalently, ^[77,77] can be expanded in powers of the source fields with expan- 
sion coefficients G,„(K[ ,K m ). The normalization factor j¥ = detC cancels the 
noninteracting vacuum diagrams, such that 5f[0,0] = in the absence of interac- 
tion. As we explain below, it will be of advantage in our case to describe the system 
not by connected Green functions but by the one-particle irreducible (1PI) vertex 
functions. Their generating functional, the effective action, is obtained from Sf by 
Legendre transformation, 

r[0,0]=^[TJ,77] + (0,T7)-(77,</.). (4) 

This functional can again be expanded in powers of the fields 0, to obtain the 
vertex functions y m {K[ ,K m ). The usual relations between and F hold, i.e., 
= 8&/8i), = 5^/577 as well as 77 = 8r/8<j), 77 = 5F/50, and 5 2 F/5050 = 
(S^/StjStj)- 1 . 



2.2 Functional flow equations 

We introduce an infrared cutoff into the bare propagator that suppresses all soft 
modes, which may be a source of divergences in perturbation theory. In bulk systems 
it is convenient to use a momentum cutoff, which suppresses momenta close to the 
Fermi surface. On the other hand, if translation invariance is spoiled by impurities 
or a particular spatial setup as for quantum dots, it is easier to use a frequency cutoff, 
excluding all Matsubara frequencies below scale A using a step function @(x): 

c^m.m^l. (5) 

'£ - hk 

This changes the microscopic model to exclude soft modes; in the limit A — > the 
original model is recovered. The easiest way to understand how the Green functions 
change with A is to derive the flow equation for the cutoff-dependent generating 
functional 

Sf A [77 , 77] = - log |./K A J Dxj/W^ fQA v y v[v ^ ] e-tolMl.uO J (6) 

with (C A )~ ! = Q A , and the normalization factor ,JV A = (detg A )~'. Taking the A 
derivative and denoting Q A = Q A , 
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d A ^ A = -d A log^Y A -— L_ / Dijrw ( W,Q A w) e^^^H'M)-^) 



where the first term, — d A \ogjV A = Ti^^C 74 ), comes from the derivative of the 
normalization factor, and the trace denotes a sum over all K. 

In the present application to quantum dots, the one-particle potential is strongly 
renormalized, and it is important to include the feedback of this renormalization 
nonperturbatively in the flow equations. This is most easily achieved in the one- 
particle irreducible (1PI) scheme, where the one-particle renormalizations are in- 
cluded in the internal propagators. Hence, we perform a Legendre transform which 
now also depends on A, 

r A [0 , 0] = <S A [t7 a , fj A ] + (0 , 77 A ) - (V , 0) , (8) 

and we find d A T A = d$ A /dA + (0 , d A 7] A ) - (d A f) A , 0) = d A ^ A after taking into 
account the derivatives also of the r\ A fields. The fundamental variables of r A are 
the fields and the 77 acquire a A dependence via the relation between and 77 
[14]. By Legendre transform of the flow Eq. (7), we obtain the 1PI flow equation 

d A r A [0,0] = d A <s A = TiQ A [c A - ( 52 Qf^ ] )~'] + (^Q A 4>) ■ (9) 

The inversion of a functional of Grassmann variables with a nonzero complex part 
is well defined as a geometric series, which involves only products of Grassmann 
variables. Therefore, we write the Hessian of F A as 

(G A )-^+f A ,[0,01, (10) 



d<j) K d<j) K i 



where the inverse full propagator (G A ) _1 = Y\ = Q A — E A according to the Dyson 
equation, and f A depends at least quadratically on the 0, fields. This decomposi- 
tion allows us to write the inverse as a geometric series, 



( -J-^L-\ =(l+G A r A r 1 G A = (l-G A f A + \G A r A } 2 ----)G / 
\d<j) K d<j) K > J V J 

which we insert into the flow Eq. (9), 

d A r A = Tr(g A [C A - G A } ) + (0 , Q A 0) + Tr(S A [f A — f A G A f A +•••])■ 



(11) 
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The first term describes the flow of the zero-point function (grand canonical poten- 
tial), while the second term expresses the change of the bare propagator with A. 
The hierarchy of flow equations is encoded in the third term which represents a loop 
with an arbitrary number of vertices f A , which contribute at least two external legs 
each, connected by full propagators G A and one single-scale propagator S A which 
has support only at frequency A , 

S A =G A Q A G A . (12) 

Given the initial conditions at A = °°, the flow equations determine r A uniquely for 
all A < oo. The flow of F A thus interpolates between the bare action r A= °° = S and 
the full solution of the problem, r A=0 = F. 



2.3 Flow equation hierarchy 

Expanding the functional flow equation for r A in powers of the source fields 0,0, 
we obtain an infinite hierarchy of coupled flow equations for the m-particle vertex 
functions y£ . This hierarchy can be represented diagrammatically and the first three 
levels are shown in Fig. 1, where each line represents the propagation of either a 
particle or a hole. For instance, the last diagram in the second line includes both 
a particle-hole and a particle-particle bubble. For our application to the SIAM it 
is convenient to distinguish these 1 , so we denote by 72(1', 2'; 1,2) the two-particle 
vertex with incoming electrons 1, 2 and outgoing electrons 1', 2'. 

The first line of Fig. 1 represents the flow equation for the self-energy E A , 

^-Z A (l',l) = -r^ e '^ + 5 A (2,2')r 2 4 (l',2';l,2), (13) 

where the labels 1,1',... are a shorthand notation for/ifi,^/ etc., and the summation 
includes Matsubara frequencies. The two-particle vertex y£ entering Eq. (13) is 
determined by the second line in Fig. 1, which in turn depends on the three -particle 
vertex etc. This infinite system of coupled differential equations can usually 
not be integrated analytically; therefore one has to resort to numerical computations 
which require to truncate the hierarchy by neglecting the flow of the higher vertices. 
This may be justified perturbatively for weak interactions since the higher vertices 
not present in the bare interaction are generated only at higher orders in the effective 
interaction. As a first step, we consider only the flow of the self-energy with all 
higher vertices remaining at their initial conditions (the bare interaction V): in our 
application, this already produces qualitatively the correct result. As a second step, 
we also include the flow of but neglect 7^, which brings us even quantitatively 
within a few percent of the known results: 



1 A derivation of the flow equation using separate particle and hole propagators from the outset 
can be found, e.g., in [20]. 
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Fig. 1 Diagrammatic representation of the flow equations for the self-energy E A , the two-particle 
vertex y£ and the three-particle vertex y£ in the 1PI formulation of the fRG. 



^-^(l',2';l ) 2) = r^£G^,5^ 4 ,f^(l',2';3 ) 4)^(3',4';l,2) (14) 



3,3' 4,4' 



-^(l',4';l,3)^(3',2';4,2) + ^(2' ! 4';l,3)^(3',l';4,2) + (3^4,3'^4') 

with the first term representing the particle-particle channel and the other two the 
particle-hole channels, respectively. The last bracket means that the two terms in the 
second line have to be repeated with the changes of variables as indicated. 

There is one technical detail which allows to simplify these equations signifi- 
cantly: the sharp frequency cutoff in Eq. (5) has the advantage that the frequency 
integrals on the right-hand side of Eqs. (13) and (14) can be performed analytically. 
The propagators contain both @(|e| — A) and <5(|e| — A) = —d J \&{\e\ — A) func- 
tions, which at first look ambiguous because the step function has a jump exactly 
where the 8 function has support. However, if one smoothes the step slightly and 
makes a change of variables from A to t = @(|e| — A) and dt = — 5(|e| — A)dA 
in the vicinity of |e| = A, the integration over the step becomes well-defined and is 
implemented by the substitution [21, 15] 
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8(x-A)f(0(x-A))->8(x-A)f f{t)dt. (15) 

Using this substitution, the product of full and single-scale propagators around the 
loop is replaced by a delta function 2 <5(|e| — A) times smooth propagators of the 
form 

G A = [Gq 1 — E A ] ~ 1 . (16) 

Note that the explicit A dependence remains only in the self-energy, such that the 
resulting propagator (and hence the flow equation) is smooth in both A and e and 
can be easily integrated numerically. 



3 Application to quantum dots 

In order to give a pedagogical example how the functional flow equations can be 
used to solve an interesting physical problem, we choose as a toy model the SIAM 
[3, 4], which is used to study, e.g., transport through a quantum dot. After integrat- 
ing out the noninteracting lead degrees of freedom (see below) this model has zero 
space dimensions, hence spin is the only quantum number, and the resulting flow 
equations are particularly simple and can even be integrated analytically if the flow 
of the two-particle vertex is neglected. Nevertheless, the physics is nontrivial, and 
the exponentially small scale (the Kondo scale) we obtain is seen neither in pertur- 
bation theory in U (which in the present example is free of infrared divergencies) 
nor in self-consistent Hartree-Fock calculations. 



3.1 The single-impurity Anderson model 

Our model consists of a single site with local Coulomb repulsion U > 0, which is 
connected to two leads I = L,R via tunnel barriers tL,R (see Fig. 2): 

H = f/(n t - ±)(nj - t)+Y J £odUo-L t i( d ° c o,i+ H - c -) > ( 17 > 

o a,l 

where c/J is the creation operator for an electron with spin a on the dot, n a = d],d a 
is the spin- <7 number operator, and c] denotes the creation operator at the end of 
the semi-infinite lead I. The leads may be modeled by tight-binding chains, Hi = 
~ t Lm=o( c l fm+\ I + h- c -) w i m hopping amplitude t. The hybridization of the dot 
with the leads broadens the levels on the dot by 17(e) = Tttfpi(e), where p;(e) is 
the local density of states at the end of lead / which we henceforth assume to be 



2 There remain additional frequency constraints if not all propagators are at the same frequency. 
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constant (wide-band limit). A local magnetic field h at the dot site splits the spin-up 
and spin-down energy levels around the gate voltage V g (oti = ±1), 

e n =V g + ah/2. (18) 

Since the leads are noninteracting, we can integrate out the degrees of freedom of 
the lead electrons in the path integral and thereby obtain a hybridization contribution 
r = Ti + ifl to the bare Green function of the dot which depends only on the sign 
of the Matsubara frequency [1], 

GoAi£) = ie-{V g + ah/2)+ir,gn{ £ ) ' (19) 

By solving the interacting many-body problem we obtain a frequency-dependent 
self-energy E a (ie) on the dot, and the full dot propagator is given by the Dyson 
equation, 

G a (ie) = [G iXa (ie)- 1 - E a (ie)]- 1 . (20) 

This self-energy could be computed by perturbation theory in the strength of the 
Coulomb interaction, e.g., in a Hartree-Fock calculation, but it turns out that this 
does not capture the physical effect that is observed in experiments and that we 
wish to describe. Instead, we will now show how the fRG can be used to compute 
with very little effort an approximation for the self-energy that leads to the correct 
physical properties of the T = conductance. 



u 




Fig. 2 A quantum dot con- 
nected to two leads (reser- 
voirs). 



In general the linear response conductance is given by the current-current cor- 
relation function (Kubo formula). At zero temperature, zero bias voltage, and for a 
single interacting site the exact conductance has the simple form [7] 

G(V g )=G ] (V g )+G l (V g ) (21) 
e 2 

G a {v g ) = -Kr Pa {o) 

in terms of the dot spectral function continued analytically to real frequencies, 

Pa(e) = --ImG CT (e + /0 + ). (22) 
n 
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3.2 Flow equation for the self-energy 



In order to implement the fRG, following Eq. (5) we introduce a cutoff in Matsubara 
frequency into the bare propagator on the dot (19), 

° ° {l£) ~ ie-(V g + ah/2)+ir Sga (e) • ( } 

This is an infrared cutoff which sets the propagator to zero for frequencies smaller 
than A (preventing these modes from being excited), but leaves the high-energy 
modes unchanged (propagating). As the cutoff scale A is gradually lowered, more 
and more low-energy degrees of freedom are included, until finally the original 
model is recovered for A — > 0. Changing the cutoff scale leads to the infinite hi- 
erarchy of flow equations for the vertex functions shown in Fig. 1. We use the flow 
equation for the self-energy (13) with the single-scale propagator substituted by 
(16). The two-particle vertex y£ is in general a complicated function of three in- 
dependent frequencies (justifying the name functional renormalization group) that 
evolves during the flow by the second line of Fig. 1. As a first approximation we 
neglect the flow of the two-particle vertex and all higher vertex functions. This can 
be justified perturbatively if the bare coupling U is small, as all the terms generated 
during the flow are of higher order in U. As a consequence, the two-particle vertex 

Y2 (iea,ie'a; iea,ie'a) = U A (24) 

remains for all A at its initial condition, U A = U A= °° = U which is the bare Coulomb 
interaction in the Hamiltonian (17). We denote a = — a. As does not depend on 
frequency in our approximation, also the self-energy does not acquire a frequency 
dependence during the flow. The flow equation for the effective level position V A = 
V g + ah/2 + L A is 

d A V A = -^ [ded(\ £ \-A)G A (ie) = - 1 ^ £ 



2k J ay 2 x e t±A ie-Vg+irsgne 

U A V A /n 



(A+rf + (v A Y 



(25) 



with initial condition V A °° = V g + oh/2 [22]. At the end of the flow, the renormal- 
ized potential V a = V A=0 determines the dot spectral function (22), 

p - (£) = ^ ( £ -y^ + r^ (26) 

which is a Lorentzian of full width 2T and height 1 jnT centered around V a - Al- 
though the true spectral function has a very different form with a very sharp Kondo 
peak at e = V a , this difference is not seen in the T = conductance (21): it probes 
only the value of the spectral function at the chemical potential e = 0, 
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G °^ = JvZTr>- (27) 

In the noninteracting case, V a — V g + oh/ 2 and the conductance is a sum of two 
Lorentzians of the applied gate voltage. If interaction is switched on, this changes 
drastically. 

In Fig. 3 the conductance G as a function of gate voltage V g for the SIAM is 
shown for different values of U /F in the upper panel, together with the occupation 
of the dot in the lower, both for the case without magnetic field, implying = 
Vj A = V A and Gf = Gj. For F «(/ the resonance exhibits a plateau [6]. In this 
region the occupation is close to 1 while it sharply rises/drops to 2/0 to the left/right 
of the plateau. Also for asymmetric barriers we reproduce the exact resonance height 
4F/Ffl/ (Fl + Fr) 2 (2e 2 /h) [1,6]. Here we focus on strong couplings U/T 3> 1. 
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Fig. 3 Upper panel: conductance as a function of gate voltage for different values of U/T and 
h = 0. Lower panel: average number of electrons on the dot. 



The solution of the differential Eq. (25) at A = is obtained in implicit form 

vJi(v)-yJo(v) Jojvg) 

vri(v)-yy (v) r (v,)' 1 ; 

with v = Vk/U, v g = V s 7t/U, 7 = T%jU , and Bessel functions /„, Y n . For \V g \ < V c 
this equation has a solution with small \V\, where v c = V c k/U is the first zero of Jq 
corresponding to V c ~ 0.77 U. For U ^>F the crossover to a solution with |V| being 
of order U (for \V g \ > V c ) is fairly sharp. Expanding both sides of Eq. (28) for small 
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|v| and |vg| gives 

V=V ? exp(-i^). (29) 

The consequent exponential pinning of the spectral weight in Eq. (22) at the chemi- 
cal potential for small \V g \ and the sharp crossover to a V of order U when \V H \ > V c 
leads to the observed resonance line shape represented by the dashed line in Fig. 4. 
For U ^> r the width of the plateau is 2V C ~ 1 .5 U, which is larger than the width 
U found with the Bethe ansatz [2, 6] corresponding to the solid line in Fig. 4. Our 
approximation furthermore slightly overestimates the sharpness of the box-shaped 
resonance. The inclusion of the renormalization of the two-particle vertex improves 
the quantitative accuracy of the results considerably (see below), while the pinning 
of the spectral function and the subsequent resonance plateau is captured already 
at the first order of the flow-equation hierarchy, even though the spectral function 
neither shows the narrow Kondo resonance nor the Hubbard bands. 




3.3 Flow equation for the fermion interaction 

With the simple approximation above, we have obtained an exponentially small 
scale for the local potential from a renormalization group flow equation. This yields 
a plateau in the conductance that qualitatively agrees already well with the exact 
Bethe ansatz solution. We will now demonstrate how one can go beyond the sim- 
plest approximations in the renormalization group flow equations to improve (sys- 
tematically for small U) the result and the agreement with the known solutions. Any 
improvement of the flow equation has to come from a more detailed parametriza- 
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tion of the two-particle vertex. The next step is to include the flow of the two-particle 
vertex y£ from the second line in Fig. 1, but still neglect the three-particle vertex y£ 
and the frequency dependence of y£ . 3 Using the vertex flow Eq. (14) with a sharp 
cutoff and the parametrization (24), we obtain [22] 



dA 



A\2 



2k 



£ G A (ie)G A (-ie) + G A (ie)G A (ie) 



£ = ±A 



2 (U A fv A V A /l 



(A+r) 2 + (v A ) 2 (A+r) 2 + {v A ) 2 



(30) 



again with initial condition U °° = U. 

A systematic improvement with respect to the previous results can be observed 
in Fig. 4, with the dotted line resulting from the inclusion of the interaction renor- 
malization. The quantitative agreement with the exact results is excellent and holds 
for U /r — 25, the largest value with available Bethe ansatz data [6]. Also for more 
complex dot geometries this extension considerably improves the agreement with 
NRG results [22]. We note that a truncated fRG scheme including the full frequency 
dependence of the two-particle vertex, and hence a frequency dependent self-energy, 
reproduces also the Kondo resonance and Hubbard bands in the spectral function 
[20]. This requires, however, a substantial computational effort. 



3.4 Effect of magnetic fields 

We next consider the case of finite magnetic fields. For h > the Kondo resonance 
in the NRG solution of the spectral function splits into two peaks with a dip at 
CO = 0, resulting in a dip of G(V g ) at V g = 0. In Figs. 5 and 6 we compare the total 
G = Gj + Gi and partial Gj conductance obtained from the above fRG truncation 
scheme including the flow of the effective interaction with NRG [23] results for 
different h expressed in units of T^ RG = 0.1 16F, where T^ RG is the width of the 
Kondo resonance at the particle-hole symmetric point V g = obtained by NRG [23]. 

The excellent agreement between NRG and fRG results provides strong evidence 
of the presence of the Kondo scale within the truncated fRG scheme. 



3.5 Determining the Kondo scale 

From the above comparison of fRG and NRG data shown in Fig. 6 for different h 
we infer the appearance of an exponentially small energy scale defining the Kondo 



3 Including the frequency dependence of y£ is an ambitious project that is beyond the scope of this 
introductory article [20]. 




Fig. 5 Gate voltage dependence of the total conductance G of a single dot with U/T = 37C 
and h/r = 0, 0.058, 0.116, 0.58 from top to bottom. In units of the V g = Kondo temperature 
7f RG /T = 0.116 these fields correspond to h = 0, 0.57^ RG , Tj? RG , and 57^ RG . Solid line: NRG 
data from [23]. Dashed lines: fRG approximation with flow of vertex. 




scale Tk(U ,V g ) by the magnetic field required to suppress the total conductance 
to one half of the unitary limit, G(U,V g ,h = Tg) = e jh. For fixed U 3> F this 
definition applies for gate voltages within the resonance plateau for h = 0. In Fig. 7 
we show Tk(U ,V g ) for different V g as a function of U on a linear-log scale. The 
curves can fitted to a function of the form 

(31) 

with Vg-dependent coefficients a, b, and c. The above form is consistent with the 
exact Kondo temperature Tg [1] that depends exponentially on a combination of U 
and the level position 
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(32) 



The prefactor of the exponential depends on the details of the model considered. To 
leading order its U and V g dependence can be neglected. For the fit to Eq. (31) we 
find b(V g ) ps 0.32 for all V g (see Eq. (29) for comparison), in good agreement with 
the exact value 7r/8 ps 0.39. The prefactor a depends only weakly on V g , and c(V g ) 
increases approximately quadratically with V g /T as shown in the inset of Fig. 7, ac- 
cording to the behavior of the exact Kondo temperature Eq. (32). We thus conclude 
that Tk{U ,Vg) can be determined from the h dependence of the conductance ob- 
tained from the fRG. Similar results are obtained from the local spin susceptibility 
[1]. 




Fig. 7 The Kondo scale Tk as a function of U for different V g : V g = (solid), V g /T — 1 (short 
dashed), V g /r = 2 (dotted), and V g /T = 4 (long dashed). Inset: The fitting parameter c as a func- 
tion of {v g /rf. 



For the single dot at T = the exact conductance, transmission phase, and dot 
occupancies are directly related by a generalized Friedel sum rule [1]: G a /(e 2 /h) = 
sin 2 (n (/iff)), and the transmission phase a a = K (n a ). As < (n a ) < 1 the argu- 
ment of sin 2 is restricted to a single period and the relation between G a , (n a ), and 
<X a is unique. In many approximation schemes the Friedel sum rule does not hold 
exactly. In contrast, within our method we map the many-body problem onto an ef- 
fective single-particle one for which the Friedel sum rule is fulfilled. For gate volt- 
ages within the h = conductance plateau the (spin independent) dot filling is 1/2 
and the (spin independent) phase is n/2. For sufficiently large U/T the crossover 
to («ff) = 1 and a<j = 7T to the left of the plateau as well as (n a ) =0 and a c = to 
the right is fairly sharp. 
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4 Summary and Conclusion 



We presented an fRG scheme developed for the study of electronic transport through 
quantum-dot systems. The resulting differential flow equations describe the effec- 
tive renormalized level position and interaction in presence of local Coulomb inter- 
actions and magnetic fields. Analytical estimates capture signatures of the Kondo 
effect, and a comparison with exact Bethe ansatz and high-precision NRG results 
shows excellent agreement up to the largest Coulomb interaction for which Bethe 
ansatz or NRG data are available. 

The presented fRG scheme constitutes a reliable and promising tool in the inves- 
tigation of correlation effects and interference phenomena in quantum-dot systems. 
The application to a single quantum dot can be directly extended to different ge- 
ometries involving more complicated setups [22] such as parallel double dots [24]. 
In [25] this method was used to investigate the long-standing phase lapse puzzle of 
the transmission phase through multi-level quantum dots. 

Furthermore, the truncated fRG scheme was sucessfully used to study one- 
dimensional fermionic lattice models with two-particle interaction and impurities 
[26, 27, 28, 29] (inhomogeneous Luttinger liquids). In addition to the universal 
low-energy asymptotics the fRG captures nonuniversal properties at higher en- 
ergy scales. Novel low-energy fixed points were found for Y-junctions of one- 
dimensional quantum wires pierced by a magnetic flux [30]. Also a single dot with 
Luttinger liquid leads [31] was investigated. For the latter the competition between 
Kondo and Luttinger liquid physics leads to a broad resonance plateau for all exper- 
imentally accessible length scales, whereas the low-energy fixed point is described 
by a sharp resonance characteristic for the Luttinger liquid behavior. 
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